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We investigate the possibility that the behavior of the rotational velocities of test particles grav- 
itating around galaxies can be explained in the framework of modified gravity models with non- 
minimal matter-geometry coupling. Generally, the dynamics of test particles around galaxies, as 
well as the corresponding mass deficit, is explained by postulating the existence of dark matter. The 
extra-terms in the gravitational field equations with geometry-matter coupling modify the equations 
of motion of test particles, and induce a supplementary gravitational interaction. Starting from the 
variational principle describing the particle motion in the presence of the non-minimal coupling, the 
expression of the tangential velocity of a test particle, moving in the vacuum on a stable circular 
orbit in a spherically symmetric geometry, is derived. The tangential velocity depends on the metric 
tensor components, as well as of the coupling function between matter and geometry. The Doppler 
velocity shifts are also obtained in terms of the coupling function. If the tangential velocity profile is 
known, the coupling term between matter and geometry can be obtained explicitly in an analytical 
form. The functional form of this function is obtained in two cases, for a constant tangential veloc- 
ity, and for an empirical velocity profile obtained from astronomical observations, respectively. All 
the physical and geometrical quantities in the modified gravity model with non-minimal coupling 
between matter and geometry can be expressed in terms of observable/measurable parameters, like 
the tangential velocity, the baryonic mass of the galaxy, and the Doppler frequency shifts. Therefore, 
these results open the possibility of directly testing the modified gravity models with non-minimal 
coupling between matter and geometry by using direct astronomical and astrophysical observations 
at the galactic or extra-galactic scale. 

PACS numbers: 04.50.+h,04.20.Cv, 95.35.+d 



I. INTRODUCTION 

The rotation curves for galaxies or galaxy clusters 
should, according to Newton's gravitation theory, show 
a Keplerian decrease with distance r of the orbital ro- 
tational speed Vt g at the rim of the luminous matter, 
vf g cx M(r)/r, where M{r) is the dynamical mass. How- 
ever, one observes instead rather flat rotation curves 
0, 0] ■ Observations show that the rotational velocities 
increase near the center of the galaxy and then remain 
nearly constant at a value of Vt g0 o ~ 200 — 300 km/s. 
This leads to a general mass profile M(r) rv t g0 o/G 
0, ■ Consequently, the mass within a distance r from 
the center of the galaxy increases linearly with r, even at 
large distances, where very little luminous matter can be 
detected. 

This behavior of the galactic rotation curves is ex- 
plained by postulating the existence of some dark (invis- 
ible) matter, distributed in a spherical halo around the 
galaxies. The dark matter is assumed to be a cold, pres- 
sureless medium. There are many possible candidates 
for dark matter, the most popular ones being the weekly 
interacting massive particles (WIMP) (for a recent re- 
view of the particle physics aspects of dark matter see 
Q). Their interaction cross section with normal bary- 
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onic matter, while extremely small, are expected to be 
non-zero, and we may expect to detect them directly. 
Scalar fields, Bose-Einstein condensates, or long range 
coherent fields coupled to gravity have also been used to 
model galactic dark matter 

However, despite more than 20 years of intense ex- 
perimental and observational effort, up to now no non- 
gravitational evidence for dark matter has ever been 
found: no direct evidence of it, and no annihilation radi- 
ation from it. 

Therefore, it seems that the possibility that Einstein's 
(and the Newtonian) gravity breaks down at the scale of 
galaxies cannot be excluded a priori. Several theoreti- 
cal models, based on a modification of Newton's law or 
of general relativity, have been proposed to explain the 
behavior of the galactic rotation curves 

A very promising way to explain the recent observa- 
tional data d, 0] on the acceleration of the Universe 
and on dark matter is to assume that at large scales 
the Einstein gravity model of general relativity breaks 
down, and a more general action describes the gravita- 
tional field. Theoretical models in which the standard 
Einstein-Hilbert action is replaced by an arbitrary func- 
tion of the Ricci scalar R, first proposed in Q , have been 
extensively investigated lately. For a review of f(R) gen- 
eralized gravity models and on their physical implications 
see @ . The possibility that the galactic dynamic of mas- 
sive test particles can be understood without the need 
for dark matter was also considered in the framework of 
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f(R) gravity models 

A generalization of the f(R) gravity theories was pro- 
posed in [15| by including in the theory an explicit cou- 
pling of an arbitrary function of the Ricci scalar R with 
the matter Lagrangian density L mat . As a result of 
the coupling the motion of the massive particles is non- 
geodesic, and an extra force, orthogonal to the four- 
velocity, arises. The connections with Modified Orbital 
Newtonian Dynamics (MOND) and the Pioneer anomaly 
were also explored, and it was suggested that the matter- 
geometry coupling may ne responsibly for the observed 
behavior of the galactic rotation curves. The model was 
extended to the case of the arbitrary couplings in both 
geometry and matter in [l6| . The implications of the 
non-minimal coupling on the stellar equilibrium were in- 
vestigated in [l7| , where constraints on the coupling were 
also obtained. An inequality which expresses a necessary 
and sufficient condition to avoid the Dolgov-Kawasaki in- 
stability for the model was derived in [18J. The relation 
between the model with geometry-matter coupling and 
ordinary scalar-tensor gravity, or scalar-tensor theories 
which include non-standard couplings between the scalar 
and matter was studied in [T|| . In the specific case where 
both the action and the coupling are linear in R the ac- 
tion leads to a theory of gravity which includes higher 
order derivatives of the matter fields without introduc- 
ing more dynamics in the gravity sector [2(| ■ The equiv- 
alence between a scalar theory and the model with the 
non-minimal coupling of the scalar curvature and matter 
was considered in [2l| . This equivalence allows for the 
calculation of the Parameterized Post-Newtonian (PPN) 
parameters /3 and 7, which may lead to a better under- 
standing of the weak- field limit of f(R) theories. The 
equations of motion of test bodies in the nonminimal cou- 
pling model by means of a multipole method were derived 
in [22|]. The energy conditions and the stability of the 
model under the Dolgov-Kawasaki criterion were studied 
in [23j]. The perturbation equation of matter on sub- 
horizon scales in models with arbitrary matter-geometry 
coupling, as well as the effective gravitational constant 
G e ff and two parameters E and 77, which along with the 
perturbation equation of the matter density are useful to 
constrain the theory from growth factor and weak lens- 
ing observations were derived in [241 ] . The age of the 
oldest star clusters and the primordial nucleosynthesis 
bounds were used in order to constrain the parameters 
of a toy model. The problem of the correct definition of 
the matter Lagrangian of the theory and of the definition 
of the energy-momentum tensor, considered in (25j and 
[III, was solved in [27j . For a review of modified f(R) 
gravity with geometry-matter coupling see [28j . 

It is the purpose of the present paper to investigate 
the possibility that the observed properties of the galac- 
tic rotation curves could be explained in the framework of 
the modified gravity theory with non-minimal coupling 
between matter and geometry, without postulating the 
existence of dark matter. As a first step in this study, 
starting from the variational formulation of the equations 



of motion, we obtain the expression of the tangential ve- 
locity of test particles in stable circular orbits. The tan- 
gential velocity is determined not only by the metric, as 
in standard general relativity, but it also depends on the 
explicit form of the coupling function between matter and 
geometry, as well as on its derivative with respect to the 
radial coordinate. In the asymptotic limit of a flat space- 
time, far away from the matter sources, the tangential ve- 
locity becomes a function of the coupling function only, 
and it does not decay to zero. Therefore the behavior 
of the neutral hydrogen gas clouds outside the galaxies, 
and their flat rotation curves, can be explained in terms 
of a non-minimal coupling between matter and geometry 
Since the observations on the galactic rotation curves are 
obtained from the Doppler frequency shifts, we general- 
ize the expression of the frequency shifts by including the 
effect of the matter-geometry coupling. Thus, at least in 
principle, the coupling function can be obtained directly 
from astronomical observations. Since the tangential ve- 
locity directly depends on the geometry matter-coupling, 
its knowledge allows the complete determination of the 
coupling function from the observational data. The form 
of the coupling function can be immediately obtained in 
the flat rotation curves region. By adopting an empiri- 
cal law for the general form of the rotation curves, and 
by adopting some reasonable assumptions on the metric, 
the coupling function can be obtained exactly for the en- 
tire galactic space-time. Therefore, all the parameters of 
the modified gravity model with linear coupling between 
matter and geometry can be either obtained directly, or 
severely constrained by astronomical observations. 

The present paper is organized as follows. The equa- 
tions of motion in modified gravity with linear coupling 
between matter and geometry, the variational principle 
for the equations of motion, as well as the matter La- 
grangian, are presented in Section II. The tangential ve- 
locity of test particles in stable circular orbits, and the 
corresponding Doppler frequency shifts are derived in 
Section III. From the study of the galactic rotation curves 
the explicit form of the matter-geometry coupling is ob- 
tained in Section IV. We discuss and conclude our results 
in Section V. In the present paper we use the Landau- 
Lifshitz 29] sign conventions and definitions, and the nat- 
ural system of units with c = I . 

II. EQUATIONS OF MOTION IN MODIFIED 
GRAVITY WITH LINEAR COUPLING 
BETWEEN MATTER AND GEOMETRY 

By assuming a linear coupling between matter and ge- 
ometry, the action for the modified theory of gravity takes 
the form 

S = J {5/1 OR) + [1 + C/ 2 (i?)] L mat \ V=g d 4 x , (1) 

where fi(R) (with i = 1,2) are arbitrary functions of 
the Ricci scalar R, and L mat is the Lagrangian density 
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corresponding to matter. The strength of the interac- 
tion between f2{R) and the matter Lagrangian L mat is 
characterized by a coupling constant (. 

By assuming that the Lagrangian density L mat of 
the matter depends only on the metric tensor compo- 
nents g^ v only, and not on its derivatives, the energy- 
momentum tensor of the matter is given by = 
L m at9iiv - 2dL mat /dg^. With the help of the held 
equations one obtains for the covariant divergence of the 
energy-momentum tensor the equation (l5l ul\ 



V"V = 2{VHn[l + Ci20R)]} 



dL r . 



(2) 



In the following we will restrict our analysis to the case 
in which the matter, assumed to be a perfect thermody- 
namic fluid, obeys a barotropic equation of state, with 
the thermodynamic pressure p being a function of the 
rest mass density of the matter (for short: matter den- 
sity) p only, so that p — p(p). In this case, the matter 
Lagrangian density becomes an arbitrary function of p, 
so that L mat — L mat (p). Then the energy- momentum 
tensor of the matter is given by 0, H?} 



dp 



dL r 



dp 



(3) 



where the four- velocity = dx^/ds satisfies the con- 
dition g^UyUv = 1. To obtain Eq. ([3]) we have im- 
posed the condition of the conservation of the matter 
current, V y {pu u ) — 0, and we have used the relation 
dp = (1/2) p {g^ — u^Uu) 6g^ v . With the use of the iden- 



tity u u W v u^ = d 2 x»/ds 2 + T[ 



from Eqs. ([2]) and 



([3|) we obtain the equation of motion of a test particle in 
the modified gravity model with linear coupling between 
matter and geometry as 



d 2 r^ 



w 



where 



/" = -V,ln [l + C/ 2 (iJ)] 



dL r , 



(p) 



dp 



(5) 

The extra- force / M , generated due to the presence of 
the coupling between matter and geometry, is perpen- 
dicular to the four-velocity, f^u^ = 0. The equation 
of motion E q. dp can be obtained from the variational 
principle (HlHl 



68,-6 1 L p ds = 6 J y/Qy/g^uWds = 0, (6) 

where S p and L p = \fQ\/gp. v u^u v are the action and the 
Lagrangian density for test particles, respectively, and 



Q = [l + tf 2 (R)} 



dLmat (p) 

dp 



(7) 



The matter Lagrangian can be expressed as [2l} 

L mat (p) = p[l + Tl(p)]- [ dp, (8) 

Jpo 

where LI (p) = dp/ p, and po is an integration con- 
stant, or a limiting pressure. The corresponding energy- 
momentum tensor of the matter is given by [271 ] 

Z" 4 " = {p [1 + * (p)] + p (p)} u»u» - p (p) <T , (9) 

respectively, where 

p(p) 



*(p)= / ^d P = n(p) 

Jpo P 



p 



(10) 



and with all the constant terms included in the definition 
of p. 



III. STABLE CIRCULAR ORBITS AND 
FREQUENCY SHIFTS IN MODIFIED GRAVITY 
WITH LINEAR COUPLING BETWEEN 
MATTER AND GEOMETRY 

The galactic rotation curves provide the most direct 
method of analyzing the gravitational field inside a spiral 
galaxy. The rotation curves are obtained by measuring 
the frequency shifts z of the 21-cm radiation emission 
from the neutral hydrogen gas clouds. Usually the as- 
tronomers report the resulting z in terms of a velocity 
field v tg fl. 

The starting point in the analysis of the motion of 
the hydrogen gas clouds in modified gravity with linear 
coupling between matter and geometry is to assume that 
gas clouds behave like test particles, moving in a static 
and spherically symmetric space-time. Next, we consider 
two observers Oe and Ooo, with four- velocities u E and 
u^, respectively. Observer Oe corresponds to the light 
emitter (i. e., to gas clouds placed at a point Pe of the 
space-time) , and Ooo represent the detector at point Poo , 
located far from the emitter, and that can be idealized 
to correspond to " spatial infinity" . 

Without loss of generality, we can assume that the gas 
clouds move in the galactic plane 8 = tt/2, so that — 

^i,f,0, (j^j , where the dot stands for derivation with 

respect to the affine parameter s. On the other hand, 
we suppose that the detector is static (i.e., Ooo's four- 
velocity is tangent to the static Killing field d/dt), and 
in the chosen coordinate system its four- velocity is = 
(4,0,0,0)^139. 

The static spherically symmetric metric outside the 
galactic baryonic mass distribution is given by 



ds 2 = e v ^dt 2 - e A Wdr 2 - r 2 (d9 2 + sin 2 Odc/i 2 ) , (11) 

where the metric coefficients are functions of the radial 
coordinate r only. The motion of a test particle in the 
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gravitational field in modified gravity with linear cou- 
pling between matter and geometry can be described by 
the Lagrangian 



L = Q 



dtY 

ds) 



Mr) 



r [ Ts) 



(12) 

where dti 2 = d9 2 + sin 2 6d4> 2 . For 9 = tt/2, dVt 2 = dcj) 2 . 
From the Lagrange equations it follows that we have two 
constants of motion, the energy E, 



E = Qe u ( r H, 
and the angular momentum I, given by 
I = Qr 2 <j>, 



(13) 



(14) 



where a dot denotes the derivative with respect to the 
affine parameter s. The condition u^u^ = 1 gives 1 = 
, from which, with the use of the 



constants of motion, we obtain 
E 2 = Q 2 e lJ+x f 2 + e v 



(15) 



This equation shows that the radial motion of the par- 
ticles is the same as that of a particle in ordinary Newto- 
nian mechanics, with velocity r, position dependent mass 
m e ff = 2Q 2 e l/+A , and energy E 2 , respectively, moving in 
the effective potential 



Veff (r) = e 



(16) 



The conditions for circular orbits dV e ff/dr = and 
r = lead to 



2 _ 1 r 3 Q {v'Q + 2Q') 



and 



E 



1 - rz//2 
, _ e-Q {rQ' + Q) 



1 - rv'/2 



(17) 



(18) 



respectively. 

The line element given by Eq. (Ill[) can be rewritten 
in terms of the spatial components of the velocity, nor- 
malized with the speed of light, measured by an inertial 
observer far from the source, as ds 2 = dt 2 (l — v 2 ) , where 



HI 



(19) 



For a stable circular orbit dr/dt = 0, and the tangential 
velocity of the test particle can be expressed as 



u tg 



-v 2 

e r 



~dt) 



(20) 



In terms of the conserved quantities the angular veloc- 
ity is given, for 9 = n/2, by 



u tg 



I^E 2 ' 



With the use of Eqs. (HTJ) and {JTHJ) we obtain 



lr{u'Q + 2Q') 
2 rQ' + Q 



(21) 



(22) 



Thus, the rotational velocity of the test body in modi- 
fied gravity with linear coupling between matter and ge- 
ometry is determined by the metric coefficient exp(^), 
and by the function Q and its derivative with respect to 
the radial coordinate r. In the standard general relativis- 
tic limit £ = 0, Q = 1, and we obtain v 2 g = rv' /2. 

The rotation curves of spiral galaxies are inferred from 
the red and blue shifts of the emitted radiation by gas 
clouds moving in circular orbits on both sides of the cen- 
tral region. The light signal travels on null geodesies with 
tangent We may restrict k* 1 to lie in the equatorial 
plane 9 — tt/2, and evaluate the frequency shift for a light 
signal emitted from Oe in circular orbit and detected by 
Ooo. The frequency shift associated to the emission and 
detection of the light signal is given by 



1 



uje_ 

Woo 



(23) 



where uii 



and the index / refers to emission 



(/ = E) or detection (/ = oo) at the corresponding 
space-time point [13, H3j- Two frequency shifts, corre- 
sponding to maximum and minimum values, are asso- 
ciated with light propagation in the same and opposite 
direction of motion of the emitter, respectively. Such 
shifts are frequency shifts of a receding or approach- 
ing gas cloud, respectively. In terms of the tetrads 
e(0) = e - v ^d/dt, e (1) = e- x i 2 d/dr, e (2) = r - 1 d/d9, 

e(x\ = (r sin 9) 1 d/d(f>, the frequency shifts take the form 



z± = 1 



where v 



T 3 



,Kc-^M]/2 



1/2 



(24) 



with uu) the compo- 



[U(i)/u (0) y 

nents of the particle's four velocity along the tetrad (i. 
e., the velocity measured by an Eulerian observer whose 
world line is tangent to the static Killing field). T = 

(l — v 2 ) is the usual Lorentz factor, and exp(j/ 00 ) 
is the value of exp[j/((r))] for r — > oo. In the case of 
circular orbits in the 9 — tt/2 plane, we obtain 



z± 



! _ c [^-v(r)l/ 2 lT y/r {v'Q + 2QQ /2 {rQ' + Q) 
y/1 - r {v'Q + 2Q>) 12 {rQ' + Q) ' 

(25) 

It is convenient to define two other quantities, z p = 
(z + — z~) /2, and za — {z+ + Z-) /2, respectively [30| . 
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In the modified gravity model with linear coupling be- 
tween matter and geometry these redshift factors are 
given by 



z D (r) 



[v. -*(r)]/3 y/r(v>Q + 2Q>)/2(rQ' + Q) 

^l~r{v'Q + 2Q') /2{rQ' + QY 

(26) 



and 



za (r) = 1 



[i^oo-^(r-)]/2 



v/l - r {v'Q + 2Q') /2 (rQ' + Q) ' 



(27) 



respectively, which can be easily connected to the obser- 
vations [301 . and Z£> satisfy the relation (z^ — l) 2 — 
z 2 D = exp [2 (z/qo — ^(r))], and thus in principle, by as- 
suming that the metric tensor component exp[f((r))] 
is known, Q and Q' can be obtained directly from the 
astronomical observations. This could provide a direct 
observational test of the galactic geometry, and, implic- 
itly, of the modified gravity models with linear coupling 
between matter and geometry. 



IV. GALACTIC ROTATION CURVES AND THE 
MATTER-GEOMETRY COUPLING 

The tangential velocity v tg of gas clouds moving like 
test particles around the center of a galaxy is not di- 
rectly measurable, but can be inferred from the red- 
shift Zoo observed at spatial infinity, for which 1 + 

Zoo = exp [(z/qc - v) /2] (1 ± v tg ) / Jl - v} g . Because of 

the non-relativistic velocities of the gas clouds, with 



v t g < (4/3) x 10 , we observe v tg 



(as the first 



part of a geometric series), with the consequence that 
the lapse function exp (v) necessarily tends at infinity to 
unity, i. e., e v « e"«"/(l-u? ff ) « e"~ ->■ 1. The ob- 
servations show that at distances large enough from the 
galactic center v tg ~ constant 

In the following we use this observational constraint 
to reconstruct the coupling term between matter and ge- 
ometry in the " dark matter" dominated region, far away 
from the baryonic matter distribution. By assuming that 
Vtg — constant, Eq. (|22[) can be written as 



<^ rQ) = 2 + - 



9L 

Q 



and can be immediately integrated to give 



Qir) - | - 



exp 



2 (1 - <) 



(28) 



(29) 



where tq is an arbitrary constant of integration. By as- 
suming that the hydrogen clouds are a pressureless dust 
(p = 0) that can be characterized by their density p 
only, the Lagrangian of the matter (gas cloud) is given 
by L mat (p) = p. Therefore, from Eqs. ([7]) and (|2T)1) we 



obtain 



caw = - 



•?,/>(■-<) 



exp 



4(l-<) 



u tg 



ln^ 



2(l-<) r 4(l-<) 



J tg 



lu- 



ll - H) r ° 



(30) 



For a general velocity profile v tg = v tg (r), the general 
solution of Eq. (j2"gj) is given by 



= l+C/ 2 (i?) = VQo exp 



1 /■ i£(r)/r - i//2 



where Qq is an arbitrary constant of integration. 



(31) 



For v 2 „ we assume the simple empirical dark halo ro- 



tational velocity law [31[ 



"0" 



J tg 



a 1 + x 



(32) 



where x — r/r optl r opt is the optical radius contain- 
ing 83% of the galactic luminosity. The parameters 
a, the ratio of the halo core radius and r opt , and the 
terminal velocity vq are functions of the galactic lumi 
nosity L. For spiral galaxies a 



1.5(£/L*) 1/5 and 



J opt 



(1 - #») (l + a 2 ), where v opt = v tg (r opt ), and 



0. = 0.72 + 0.441og 10 (L/L*), with = 1O 1O - 4 L . 

By assuming that the coupling between the neutral hy- 
drogen clouds and the geometry is small, C/2 (R)L ma t « 

1, and consequently does not significantly modify the 
galactic geometry, the vacuum metric outside the bary- 
onic mass distribution with mass Mb, corresponding to 
Lmat ~ 0, is given by the static spherically symmet- 
ric solution of the f(R) modified gravity. By assuming, 
for simplicity, that the galactic metric is given by the 
Schwarzschild metric (which is also a solution of the vac- 
uum field equations of the f(R) modified gravity (HJ), 
written as 



1 



2-Rp 

: 

X 



(33) 



where Rq = GMb /r opt , from Eq. (|3"Tj) we obtain 



I + C/2OR) = exp 



a x arctanh — -x 



/ 2R Q y [{l-vl) x 2 + a^ 
V x J Qn 1/2 x x ^ 4 



Qo x 



where 



qRqVq 



2^vlW+A{l-vl)Riy 



(34) 



(35) 
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ml 



4 [a 2 + 4 (1 



and 



vt[a* + Z(l-vl)Rl] 
±{l-vl)[a?+A{l-vl)Rl] 



(36) 



(37) 



respectively. Thus the geometric part of the coupling 
between matter and geometry can be completely recon- 
structed from the observational data on the galactic ro- 
tation curves. 



V. DISCUSSIONS AND FINAL REMARKS 

The galactic rotation curves and the mass distribution 
in clusters of galaxies continue to pose a challenge to 
present day physics. One would like to have a better un- 
derstanding of some of the intriguing phenomena associ- 
ated with them, like their universality, and the very good 
correlation between the amount of dark matter and the 
luminous matter in the galaxy. To explain these observa- 
tions, the most commonly considered models are based 
on particle physics in the framework of Newtonian grav- 
ity, or of some extensions of general relativity. 

In the present paper we have considered, and further 
developed, an alternative view to the dark matter prob- 
lem [15| , namely, the possibility that the galactic rotation 
curves and the mass discrepancy in galaxies and clusters 
of galaxies can naturally be explained in gravitational 
models in which there is a non-minimal coupling between 
matter and geometry. The extra-terms in the gravita- 
tional field equations modify the equations of motion of 
test particles, and induce a supplementary gravitational 
interaction, which can account for the observed behav- 
ior of the galactic rotation curves. As one can see from 
Eq. (f2"2"j). in the limit of large r, when v' — > (in the case 
of the Schwarzschild metric v' w 2Mb/t 2 ), the tangential 
velocity of test particles at infinity is given by 



J tgoo 



Q/Q" 



(38) 



which, due to the presence of the matter-gravity coupling, 
does not decay to zero at large distances from the galactic 
center, a behavior which is perfectly consistent with the 
observational data [![, and is usually attributed to the 
existence of the dark matter. 

By using the simple observational fact of the constancy 
of the galactic rotation curves, the matter-geometry cou- 
pling function can be completely reconstructed, without 
any supplementary assumption. 

If, for simplicity, we consider that the metric in the 
vacuum outside the galaxy can be approximated by the 
Schwarzschild metric, with exp (v) = 1 — 2GMs/r, where 
Mb is the mass of the baryonic matter of the galaxy, 
then, in the limit of large r, we have v — > 0. Therefore 



from Eq. (|30j) we obtain 



C lim h (R) 



J tg 



(l-< 



■In- 



fo 



(39) 



If the galactic rotation velocity profiles and the galactic 
metric are known, the coupling function can be recon- 
structed exactly over the entire mass distribution of the 
galaxy. 

One can formally associate an approximate " dark mat- 
ter" mass profile Mdm (r) to the tangential velocity pro- 
file, which is determined by the non-minimal coupling 
between matter and geometry, and is given by 



M DM {r) 



1 r 2 (v'Q + 2Q') 
2G rQ' + Q 



(40) 



The corresponding "dark matter" density profile 
Pdm ( r ) can be obtained as 



Pdm (r) 



1 dM 



1 



Airr 2 dr AitG 
v'Q + 2Q' v"Q + vQ' + 2Q" 



r (rQ' + Q) 2 (rQ' + Q) 

VQ + 2Q') {rQ" + 2Q')~ 
2 (rQ' + Qf 



(41) 



In order to observationally constrain Mdm and pdm 
we assume that each galaxy consists of a single, pressure- 
supported stellar population that is in dynamic equilib- 
rium and traces an underlying gravitational potential re- 
sulting from the non-minimal coupling between matter 
and geometry. Further assuming spherical symmetry, the 
equivalent mass profile induced by the geometry - matter 
coupling (the mass profile of the " dark matter" halo) re- 
lates to the moments of the stellar distribution function 
via the Jeans equation 0| 



d 
dr 



k 2 >] 



2p a (r)P _ G Ps M DM (r) 



(42) 



where p s (r), (v 2 ), and /3(r) = 1 — (w 2 ) / (v 2 ) describe 
the three-dimensional density, radial velocity dispersion, 
and orbital anisotropy of the stellar component, where 
(rig) is the tangential velocity dispersion. By assuming 
that the anisotropy is a constant, the Jeans equation has 
the solution 133] 



Ps (v 2 ) = Gr- 2 ? / ^-Vp. (s) M DM (a) da. (43) 



With the use of Eq. (I4U1) we obtain for the stellar velocity 
dispersion the equation 

p (v 2 )* 1 -^ r^p (a) ^WW+^W j, 
*W~2 f I S Ps[S) S Q'( S )+Q(s) dS - 

(44) 
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Projecting along the line of sight, the "dark matter" 
mass profile relates to two observable profiles, the pro- 
jected stellar density I(R), and to the stellar velocity 
dispersion a p {R) according to the relation |2| 



i(R) J R 



1-/3 



R u 



Ps{vl)' 

vr 



R 2 



-.dr. 



(45) 



Given a projected stellar density model I(R), one re- 
covers the three-dimensional stellar density from p s (r) = 

-(l/ir)f™{dI/dR)(R 2 -r 2 y 1/2 dR @. Therefore, 
once the stellar density profile I(R), the stellar veloc- 
ity dispersion (v 2 }, and the galactic metric are known, 
with the use of the integral equation Eq. (|44l) one can 
obtain the explicit form of the geometry-matter cou- 
pling function Q, and the equivalent mass profile in- 
duced by the non-minimal coupling between matter 
and geometry. The simplest analytic projected den- 
sity profile is the Plummer profile |2j, given by I(R) = 

C ( nr lal f) (l + R2 / r lal /) > Where £ is the t0tal lu " 

minosity, and r^ai / is the projected half-light radius (the 
radius of the cylinder that encloses half of the total lu- 
minosity). 

An important physical requirement for the circular 
orbits of the test particle around galaxies is that they 
must be stable. Let rg be a circular orbit and con- 
sider a perturbation of it of the form r = ro + 5, where 
5 << ro (34[. Taking expansions of V e ff (r), exp (y + A) 
and Q 2 (r) about r = r , it follows from Eq. (fT5)) that 



(46) 



The condition for stability of the simple circular orbits 
requires V"tf {ro) > (34|. This gives for the coupling 
function Q the constraint 



d 2 Q 



dr 2 



< 



Q' 



2l 2 A 6l 2 

r i T + 2QQ'\ - — 



(47) 



a condition that must be satisfied at any point ro of the 
galactic space-time. 

From observational as well as from theoretical point 
of view an important problem is to estimate an upper 
bound for the cutoff of the constancy of the tangential 
velocities. If in the large r limit the coupling function 
satisfies the condition Q'/Q —> 0, then Q/Q' —> oo, and 
in this limit the tangential velocity decays to zero. If 
the exact functional form of Q is known, the value of 
r at which the rotational velocity becomes zero can be 
accurately estimated. 

In the present model all the relevant physical quanti- 
ties, including the "dark mass" associated to the cou- 
pling, and which plays the role of the dark matter, 
its corresponding density profile, as well as the matter- 
geometry coupling function, are expressed in terms of 
observable parameters - the tangential velocity, the bary- 
onic (luminous) mass, and the Doppler frequency shifts 
of test particles moving around the galaxy. Therefore, 
this opens the possibility of directly testing the modified 
gravity models with non-minimal coupling between mat- 
ter and geometry by using direct astronomical and as- 
trophysical observations at the galactic or extra-galactic 
scale. In this paper we have provided some basic theo- 
retical tools necessary for the in depth comparison of the 
predictions of this model and of the observational results. 
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